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Using  four  common  catalysts  the  existence  and  magnitude  of  synergistic 
effects  upon  burn  rate  are  studied  with  HTPB-AP  composite  solid  propellants. 
The  effect  under  study  is  one  whereby  at  a fixed  total  catalyst  loading  the 
mixture  of  two  catalysts  would  be  more  effective  in  burn  rate  augmentation 
than  a single  catalyst.  At  2000  psia  all  possible  combinations  of  catalysts 
showed  synergistic  effects  but  the  effects  are  weak  and  within  the  reproduci- 
bility range  of  the  experiments.  Ferric  oxide  and  iron  blue  showed  synergistic 
effects  only  at  300  and  2000  psia  over  the  pressure  range  30C-2000  psia. 

Analysis  of  the  effect  of  b*  n4er-ojddiser  reactions  upon  the  surface 
shape  of  a two-dimensional  sandwich  has  shown  that  a)  the  faster  are  those 
reactions  the  more  deeply  will  the  sandwich  obtain  a V-aiiape  and  b)  these 
reactions  are  generally  fast  chough  to  indeed  affect  the  surface  shape. 
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INTRODUCTION 


The  current  investigation  represents  a continuation  of  the  work  initiated 
at  the  Georgia  Institute  of  Technology,  which  has  been  presented  in  References 
1 and  2.  In  these  investigations  two-dimensional  composite  solid  propellant 
sandwiches  of  polycrystalline  ammonium  perchlorate  (AP)  oxidizer  and  hydroxyl 
terminated  polybutadine  (HTPB)  binder  were  used  to  investigate  various  modes  of 
catalytic  behavior.  Four  proven  burn  rate  modifiers  which  behave  as  catalysts 
have  been  used  throughout  these  ongoing  investigations.  Hiey  are  Harshaw  Catalyst 
Cu  0202  P,  ferric  oxide,  iron  blue  and  ferrocene. 

During  tiie  first  phase  of  this  investigation^  the  catalysts  were  loaded 

in*o  either  the  oxidizer  or  binder  or  restricted  to  the  binder-oxidizer  interface 

of  two-dimensional  sandwiches.  Vertical  sandwich  burn  rates  and  oxidizer  normal 

regression  rates  were  obtained  using  cinephotomacrography  techniques.  Details  of 

the  surface  structure  were  obtained  by  using  the  scanning  electron  microscope  to 

observe  partially  burned  samples.  It  was  determined  that  if  one  effective  catalyst 

was  added  to  the  oxidizer  and  another  to  the  binder,  there  was  a net  increase  of 

burn  rate  over  that  of  the  sum  of  the  independent  actions  of  the  catalysts,  lids 

combined  effect  was  denoted  as  a positive  synergistic  effect  on  sample  burn  rate. 

The  existence  of  tlds  synergistic  effect  for  the  four  catalysts  was  systematically 

(2) 

investigated  using  the  two-dimensional  sandwiches.  The  Harshaw  catalyst  Cu  0202  P 
ferric  oxide  system,  Pig.  1,  and  iron  blue-ferric  oxide  system,  Fig*  2,  were  found 
to  exhibit  a positive  synergistic  effect  over  the  entire  pressure  range,  6t>0  to 
2000  psia  and  the  tsaximun  synergistic  effect  at  600  psia,  respectively.  The  burn 
rate  ratio  on  Figs.  1 and  2 is  the  ratio  of  the  burn  rate  of  the  catalysed  sam- 
ple to  the  burn  rate  of  a pure  AF-HTPB  sample. 
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Synergistic  Effect  on  Sandwich  Vertical 
Burn  Bate  for  Ferric  Oxide  and  Iron  Blue 


These  two  systems  were  chosen  for  farther  testing  in  cast  propellant  con- 
figuraLions « The  initial  samples  were  prepared  with  the  same  volumetric  loading 
as  in  the  sandwiches . The  strands  were  prepared  from  a common  lot  of  uncured 
composite  propellant  with  an  82.7%  solids  to  17*3%  binder  loading.  The  oxidizer 
particle  distribution  was: 

180  am  < 70%  < 212  am 
37  am  < 30%  < 45  am  . 

The  sample  burn  rate  for  the  uncatalyzcd  propellant  agreed  with  the  data  for  pro- 
pellant if  78  in  the  Princeton  University  test  series  ^ . The  catalyst  was  added 
to  the  uncured  propellant  at  a weight  percent  of  2.41.  Tests  of  the  Harshaw  cat- 
alyst Cu  0202  P - ferric  oxide  and  iron  blue-ferric  oxide  systems  showed  no  pos- 

(2) 

it'1" va  synergistic  effect  at  600  psiav  ' . 

These  cost  propellant  samples  were  initially  prepared  to  contain  on  average 
the  same  volumetric  loading  of  catalyst  as  the  sandwiches.  But  for  the  cast  pro- 
pellants all  of  the  catalyst  is  suspended  in  the  binder  matrix,  which  constitutes 
17V f>  of  the  sample.  Since  2.4y%  of  catalyst  was  added  to  the  cast  propellant  saa- 
' ‘s  the  hinder  matrix  contained  l4.2v%  of  the  catalyst.  This  binder  matrix  for  the 
sandwiches  contained  4.4%  of  catalyst.  This  increase  of  catalyst  in  the  binder 
affects©  the  eyre  and  the  strength  of  the  binder.  This  high  catalyst  loading  in 
the  binder  may  be  above  the  amount  required  for  maximum  burn  rate  augmentation. 

In  order  to  investigate  tills  possibility  further,  a new  set  of  samples  were  pre- 
pared using  a reduced  total  catalyst  loading  of  1%  catalyst  by  weight.  Tills  re- 
duced the  effective  binder  loading  to  §.<?%. 

TSte  current  study  has  continued  the  investigation  of  synergistic  effects  on 
sample  burn  rate  in  cast  composite  propellants  with  lv%  of  catalyst.  An  experimental 
procedure  has  been  developed  to  produce  a cast  propellant  that  can  be  used  as  a con- 
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vonient  tost  vehicle  for  catalytic  effects.  An  investigation  of  catalyst  loading  j 

in  the  oxidizer  of  a cast  propellant  has  been  performed.  Burn  rates  of  those  aara-  j 

pies  have  been  obtained  at  600  and  2000  psia.  ! 

1 

FUSE  WIRE  TESTS 

IWo  separate  batches  of  AP  and  IITPB  were  prepared  for  the  initial  tests  of 
the  addition  of  lw$  catalyst  to  the  uncured  propellant.  The  ratio  of  82. 7%  solids 
(AP)  to  17.3$  binder  ( IITPB)  was  maintained,  consistent  with  earlier  tests v i 

The  oxidizer  particle  distribution  was  maintained  at  : 

180  MJU  < 70fj  < 212  pm 

37  urn  < 30$  < 45  om  . 

Freon  T F was  mixed  with  the  first  batch  of  propellant,  This  decreased  the  vis- 
cosity of  the  mixture  and  allowed  the  propellant  to  be  packed  much  easier  in  the 
teflon  strand  molds.  Trie  Freon  T F could  be  removed  by  vaporisation  at  room  tem- 
perature or  by  vacuum  in  the  curing  oven.  The  samples  were  not  sufficiently  fluid 
to  fill  the  voids  as  the  freon  vaporised.  Tills  yielded  porous  cast  propellants. 

All  of  the  oxidizer  particles  were  prefixed  before  the  binder  and  catalyst  were 
added  to  the  second  batch  of  propellant.  This  propellant  was  very  difficult  to 
mix  and  pack  in  the  teflon  colds.  These  samples  resembled  the  first  batch  of 
samples.  Uiey  were  porous  with  a grainy  appearance. 

Thirteen  fuse  wire  tests  were  conducted  using  samples  froa  these  two  batches 
of  propellants.  Several  variations  of  sample  preparations  were  used,  Different- 
codes  of  inhibition  were  tried  on  the  sides  of  the  samples.  Coatings  of  FVC, 
vaccus  grease  and  epoxy  were  used  along  with  water  leachisig  of  surface  A?  and 
increased  flushing  nitrogen  flows.  The  results  for  all  samples  were  the  same. 

All  fuse  wires  appeared  to  break  simultaneously.  Both  batches  of  propellant 
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behaved  similarly,  This  had  not  been  encountered  with  the  2.4lw$  catalytic  samples 
of  Reference  2.  There  had  been  fluctuations  of  burn  rates  but  not  consistently 
high  burn  rates  as  indicated  by  these  samples. 

High  speed  motion  pictures  were  taken  of  four  samples  from  these  two  batches. 
Based  on  a frame  rate  of  1600  frames  per  second,  the  pure  AP-HTPB  sample  and  the 
1$  Harsh&w  catalyst  Cu  0232  P sample  exploded  and  burned  in  .015  seconds  (24  frames) 
The  1 $ ferric  oxide  sample  was  completely  consumed  in  .009  seconds  (14  frames) . 

The  1/u  ferric  oxide-iron  blue  sample  lasted  for  .022  seconds  (36  frames) . 

It  is  believed  that  during  manufacturing  too  narrow  an  AP  particle  size  dis- 
tribution was  being  eznployed,  preventing  binder  wetting,  and  good  packing.  It  was 
decided  to  broaden  the  particle  size  distribution.  Based  on  the  length  of  time 
fciiat  was  necessary  to  prepare  fuse  wire  samples  and  the  possible  areas  of  un- 
certainty when  the  sample  burns  rapidly,  it  was  decided  to  return  to  the  cine- 
photomnerography  for  determination  of  sample  burn  rate  and  uniformity  of  burn 
rate.  Poor  quality  propellant  samples  can  be  determined  immediately  from  viewing 
the  motion  pictures. 


ClBRPHOTOKftCROCRAFirf 


The  samples  for  these  tests  were  prepared  from  coason  lots  of  uncured  pro- 
pellant. The  solids  to  binder  ratio  was  maintained  at  82*7$  to  1?.3$.  The  bi- 
m4&l  AP  particle  distribution  was  broadened  to: 

lasa  < 35$  212  ym 

125  uffl  < 3%  < USo  ua 

k’j>  ym  < 10$  < 63  Ji«3 

37  m < 10$  < 4$  y»  j :p$>  < 63  ma 
10£  < 37  ua. 


yea 

3 


125  y»  < 70%  < 212  ys 


»- 
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The  oxidizer  particle  si2e  distribution  was  prepared  from  ulta -pure 
runmonium  perchlorate  supplied  by  the  Naval  Weapons  Center  at  China  Lake,  Calif. 

It  was  ground  in  a 0.3  gallon  porcelain  mill  jar  with  an  assortment  of  1.3  and 
P.1}  cm  (*"  and  1")  diameter  balls.  A series  of  seven  20.3  cm  (8")  diameter  U.3.A. 
standard  sieves  was  used  to  se]>arate  the  ground  AP.  A motor  driven  sieve  shaker 
was  used  to  vibrate  the  sieves.  The  ground  and  sieved  AP  was  stored  in  a sealed 
container  and  dried  for  24  horn’s  at  70°  C before  mixing  with  the  binder. 

The  propellant  was  prepared  by  mixing  all  of  the  ingredients  of  the  hydroxyl 
terminated  polybutadiene  with  approximately  one  half  of  the  largest  size  AP  par- 
ticles. Tills  was  mixed  throughly  and  all  of  the  AP  particles  are  vetted  with 
binder.  This  process  was  repeated  with  all  AP  particle  sizes  adding  the  finest 
AP  last,  litis  propellant  was  then  divided  for  catalyst  addition.  A mi m. mum  of 
10  grams  of  propellant  was  needed  for  each  catalyzed  sample.  The  proper  amount 
of  catalyst  was  added  to  the  propellant  and  the  mixing  process  repeated.  This 
propellant  was  cast  in  22  mm  by  22  mm  embedding  colds  and  cured  for  seven  days 
at  60 °C. 

The  physical  properties  and  appearance  of  the  cured  propellant  was  improved 
fey  the  increased  mixing  time  of  the  propellant  and  the  broadened  AP  particle 
distribution.  All  of  the  oxidizer  appeared  to  be  wetted  with  binder.  The  cured 
propellant  was  cut  into  convenient  sixes  ('  mrs  x Sfts?  x Igmrs)  for  testing  in  the 
high  pressure  combustion  apparatus  designed  by  Varney  ^^5  Jones  . 

These  samples  were  similar  In  si so  to  the  two-disaensionel  sandwiches  tested 
previously*  [Tils  allowed  the  sane  cinepho tecac nagrapfcl c tedmiques  to  be  used  to 
observe  the  burning  samples.  The  sane  nitrogen  flows  that  were  determined  experi- 


mentally  for  the  sandwiches  were  used  for  smoke  control.  The  latent  image  mag- 
nification was  reduced  to  slightly  loss  than  1:1  (image:  actual),  'the  motion 
picture  speed  was  maintained  at  1600  frames  per  second,  The  sample  was  oriented 
on  the  sample  holder  to  allow  observation  of  two  sides  of  the  sample.  iMs  indi- 
cated the  uniformity  of  the  burn  arid  any  effect  of  the  light  source  on  burn  rate. 

The  sample  burn  rates  and  their  uniformity  were  determined  by  recording  sam- 
ple profiles  at  100  frame  intervals  by  tracing  the  projected  burning  surface,  These 
could  be  used  to  determine  a burn  rate  every  .0625  seconds. 

The  same  combination  of  catalysts,  Harshav  catalyst  Cu  0202  P - ferric  oxide 
and  iron  blue  - ferric  oxide,  that  wer . four.' to  exhibit  positive  synergistic  ef- 
fects on  sandwich  vertical  burn  rate,  vere  tested  at  Iw $ addition  to  cast  propel- 
lants. These  two  combinations  of  catalyst  were  determined  by  the  two-dimensional 

(2) 

screening  experiment  and  wore  tested  at  2.Ul catalyst  addition  at  600  psiav  ‘ . 


Initially  the  six  samples  were  tested  at  6QQ  psia.  Ko  positibe  synergism  was 
indicated,  so  a second  set  of  samples  were  tested  at  2000  psia.  The  results  of 
both  of  these  tests  are  shown  in  table  I. 

■mis  1 

Sura  Kates  for  Cast  Composite  Propellant  (il«22-?k) 


Catalyst 


Burn  Kate  (in/sec) 

£ga  psist  20OD  psi 


r 

lior.e 

.3^ 

Harshav  Catalyst  Cu  0202  P (CC)  . 

Perris  Oxide  Cl©} 

M 

Iren  Hue  (IP-) 

.$k 

oc  & 10 

M 

10  St  IB 

-5’*» 

* * I possible  1 

I.59  Js^icnris^  |.?8 


A "positive”  synergistic  effect  on  sample  burn  rate  is  obtained  when  the 
burn  rate  of  the  sample  with  two  effective  catalysts  present  exceeds  the  burn 
rates  of  both  of  the  corresponding  samples  with  only  one  catalyst  present.  Hie 
total  catalyst  loading  is  maintained  constant  for  all  three  samples.  A "possible" 
synergistic  effect  on  sample  burn  rate  is  obtained  when  the  burn  rate  of  the 
sample  with  two  catalysts  present  is  greater  than  the  average  of  the  burn  rates 
for  the  two  samples  with  only  one  catalyst  present.  Consequently,  a "positive" 
synergism  is  also  a "possible"  synergism.  A "negative"  synergistic  effect  on 
sample  burn  rate  is  obtained  when  the  burn  rate  of  the  sample  with  two  catalysts 
present  is  less  than  the  burn  rates  of  both  of  the  corresponding  samples  with 
only  one  catalyst  present. 

The  Harshaw  catalyst  Cu  0202  P - ferric  oxide  system  did  not  exhibit  a 
positive  synergism  at  either  pressure,  but  there  was  a possible  synergism  in- 
dicated at  2000  psi.  There  was  a positive  synergistic  effect  on  sample  burn  rate 
for  the  ferric  oxide  - iron  blue  system  at  2000  psi.  Because  of  this  the  ferric 
oxide  - iron  blue  samples  and  their  companion  samples  were  tested  over  the  pres- 
sure range  300  to  2000  psia.  The  results  are  shown  in  Figure  3.  There  was  a pos- 
itive synergistic  effect  at  300  and  2000  psia,  a possible  synergistic  effect  at 
600  psia,  and  a negative  effect  at  1000  and  1500  psia.  The  positive  synergistic 
-effect  on  sample  burn  rate  at  300  and  2000  psia  along  with  a review  of  the  favor- 
able catalytic  effect  of  ferrocene  when  it  was  restrained  to  the  binder-oxidiscr 
interface  ^ led  to  a new  evaluation  of  possible  synergistic  effects  for  all  com- 
binations of  catalysts.  Five  new  catalyst  configurations  were  prepared  to  complete 
the  test  matrix.  This  included  a sample  with  1$  ferrocene  added,  three  samples  with 
ifo  of  ferrocene  and  of  the  other  three  catalysts  and  finally  a §$  Harshaw  catalyst 
Cu  0202  P and  | $ ircr.  blue  sample.  The  results  are  shown  in  Ihble  II  for  2000  psia. 
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c Effect  for  Ferric  Qxid 
st  Composite  Propellants 


TABLE  II 

Burn  Kates  for  Cast  Composite  Propellant  (5-8-75) 


Catalys  t 

2000  psi 

Burn  Rate  (in/sec) 

r/r 
' pure 

None 

.80 

1 

Ferrocene  (f) 

.75 

.94 

CC  & IB 

1.13 

1.41 

CC  & F 

.95 

1.19 

10  & F 

.89 

1.11 

IB  & F 

•87 

1.09 

Initially  these  turn  rates  were  ratloed  with  the  pui-o  cast  propellant  burn  rate 
of  Table  1 (11-22-74) . These  burn  rate  ratios  are  shown  in  the  first  column  of 
Table  III,  A positive  synergistic  effect  on  sample  burn  rate  was  indicated  for 
all  new  catalyst  combinations.  This  meant  that  only  the  Harshaw  catalyst  Cu  0202  P - 
ferric  oxide  system  did  not  yield  a positive  synergistic  effect  at  2000  psi.  Six 
separate  determinations  of  sample  burn  rate  yielded  a different  average  burn  rate 
of  the  pure  cast  composite  pi’opellant  (5-8-75)  ♦ When  this  new  burn  rate,  r^^  = 

.30  in/sec  was  used  to  ratio  the  burn  rates  of  samples  prepared  from  the  propel- 
lant (5-8-75) , the  second  column  of  burn  rate  ratios  was  obtained  in  Table  III. 

The  initial  indication  of  five  positive  syr.ergisma  was  reduced  to  positive  syner- 
gisms  for  thre*  systems. 


Due  to  the  apparent  large  difference  of  burn  rate  for  the  pure  cast  composite 
propellant  samples  a third  batch  of  propellant  was  prepared.  All  combinations  of 
the  four  catalysts  and  their  comparison  samples  were  prepared.  U;»  burn  rates  and 
their  ratios  are  shown  in  Tbble  IV  for  2000  psia.  There  was  a positive  synergistic 
effect  indicated  for  the  iron  blue  - ferrocene  system.  All  indications  of  syner- 
gistic effects  were  different  from  t!ic  ether  two  baschc : of  prcpullq.t.  There 
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TABLE  III 

Burn  Rate  Ratios  for  Cast  Composite 

Propellants  (11-22-74)  and  (5-8-75) 
at  2000  psia 


Catalyst 


i 

HarsLav  Catalyt  ; Cu  02Q2P(TC  1.12 
Ferric  Oxide  (10)  1.26 

Iron  Blue  (13)  1.16 

Ferrocene  (f) 


CC  & 10 
CC  & IB 
CC  h F 
10  & IB 


10  & F 

IB  & P 


r/r 

' pure 

rpure  = .66  in/sec 

1 

1.12 

1.26 

1.18 

1.14 

j Possible 

1 Synergism 

' pure  pure  = .66  in/sec 
or  ,80  in/sec 


appears  to  be  a spread  of  up  to  .26  in/sec  in  the  burn  rate  of  the  pure  cast 
samples . This  variation  is  as  large  as  the  burn  rate  augmentation  obtained 
for  most  of  the  catalyzed  samples  both  with  single  and  double  catalysts  present. 
It  is  concluded  that  reproductibility  with  these  small  mixes  prevents  a conclu- 
sive demonstration  of  synergistic  effects. 

TABLE  IV 

Burn  Rates  for  Cast  Composite 
Propellant  (8-20-75) 


Catalyst 

r (in;  v 
' /sec) 

r/r 
' pure 

None 

•56 

1 

Harshaw  Catalyst  Cu  0202  P (cc) 

•78 

1.39 

Ferric  Oxide  (10) 

• 

GO 

1.52 

Iron  Blue  (IB) 

•78 

1.39 

Ferrocene 

•65 

1.16 

Possible 

Synergism 

CC  & 10 

•77 

1.38 

No 

CC  & IB 

.78 

1.39 

• 

CC  & F 

.62 

1.11 

No 

10  S:  IB 

.81 

1.45 

No 

10  & F 

.67 

1.20 

No 

IB  & F 

.92 

1.64 

Yes 

CATALYST  LOADED  INTO  AP 

liars!  taw  catalyst  Cu  0202  P van  selected  for  tests  of  loading  in  the  oxidizer 
polycrystalline  structure  because  of  its  catalysis  of  the  AP  deflagration  process 
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The  oxidizer  particles  for  the  investigation  of  catalyst  loaded  into  the 
oxidizer  were  prepared  from  special  polycrystalline  oxidizer  disks.  Ultra-pure 
ammonium  perchlorate  less  than  37  pmi  in  diameter  was  mixed  with  lw$  of  catalyst 

(k) 

and  placed  in  the  mold  assembly  designed  by  Varney  * ' . Approximately  5 grams 
of  this  fine  oxidizer  mixture  was  pressed  at  one  time,  The  polycrystalline  disk 
was  subjected  to  30,500  psia  for  one  hour.  These  disks  were  reground  by  hand  in 
a motar  and  pestle  and  then  sieved.  All  oxidizer  particles  less  than  37  u™  in 
diameter  were  discarded.  10$,  of  all  remaining  sizes  was  retained  for  a catalyst 
loading  determination.  The  remaining  oxidizer  with  catalyst  was  cast  into  a 
solid  propellant.  Die  oxidizer  particle  distribution  was: 

125  um  < 43.7$  < 212  pm 
63  M®  < 50.8$  < 125  M*m 
37  um  < 5.5$  < 63  Mm 

Comparison  samples  of  pure  AP-HTPB  and  lw $ of  catalyst  added  to  the  AP-HTPB 
mixture  were  prepared  with  the  same  oxidizer  particle  distribution. 


These  three  samples  were  tested  at  600  and  2000  psia.  The  results  are  shown 
in  Table  V. 


TABLE  V 

i lor shaw  Catalyst  Cu  0202  P Loaded  in  AP 


Burn  Rate 

(in/sec) 

Propellant 

600  psia 

2000  psia 

Pure  AP-HTPB 

.82 

1$  CC  in  AP-Pure  HTPB 

•54 

.86 

1$  CC  in  HTPB- Pure  AT 

.69 

1.44 

Die  maximum  burn  rate  was  obtained  for  the  catalyst  added  to  the  AP-iiTPB  mixture 
This  is  essentially  a binder- loaded  cast  composite  propellant. 


All  three  of  the  special  distribution  propellants  were  difficult  to  mix.  Die 
cured  propellants  were  grainy  and  crumbly  in  sections.  Two  separate  batches 
of  the  pure  propellant  were  prepared  and  the  burn  rate  varied  as  .82  + .05  in/sec. 

A procedure  for  deterraing  the  amount  of  catalyst  remaining  in  the  polycrys- 
talline oxidizer  structure  was  developed.  The  oxidizer  mixture  was  placed  in  a 
60  ml  vacuum  funnel  with  a built-in  40  ram  diameter  ultra- fine  fritted  glass  fil- 
ter disk.  Die  ammonium  perchlorate  was  dissolved  by  repeated  washings  with  dis- 
tilled water.  Die  washing  was  continued  for  three  cycles  after  the  liquid  showed 
no  perchlorate  ions  when  tested  with  a methylene  blue  indicator  solution.  Diis 
solution  turns  from  dark  blue  to  violet  in  the  presence  of  the  perchlorate  ions. 

The  material  remaining  in  the  filter  wes  dried  and  weighed.  A catalyst  deter- 
mination of  the  10$  of  all  size  distributions  uaad  in  the  propellant  was  perform- 
ed along  with  the  oxidizer  with  a diameter  lean  than  37  um.  Dais  allowed  a catalyst 
balance  to  be  performed. 

Hi U ally  1 ,w$  of  Harshaw  catalyst  was  added  to  the  AP.  The  measurement  in- 
dicated a loading  of  1.6  ± .3 % by  weight  in  the  AP  which  went  into  the  propellant. 
A catalyst  balance  was  made  by  determining  the  amount  of  catalyst  in  the  sieved 
polyeivsthLU.no  oxidizer  loss  than  37  wu  in  diameter.  Diis  balance  indicated  a 
possible  oxidizer  loading  of  .7  ± .2$  by  weight.  The  difference  between  these  two 
indicated  loadings  is  unresolved. 


ANALYSIS 


ho  work  of  Kef.  (2)  was  extended  to  include  the  effect  of  binder-oxidizer 
reactions  on  the  surface  shape  attained  by  a two-dimensional  sandwich,  assuming 
there  are  no  melt  flows . Die  work  is  presented  in  Appendix  A.  It  is  shown  that. 
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indeed,  fast  binder-oxidizer  reactions  will  force  a sandwich  to  go  to  a V-shaped  f 

structure,  as  has  been  suspected  from  experimental  observation  and  physical  rea-  | 

I 

soning.  More  importantly,  it  appears  that  the  oxidizer-binder  kinetics  are  fast  1 

I 

f 

enough  to  influence  the  sandwich  shape  at  usual  deflagration  pressures  and  that  j 

the  AP  flame  is  not  the  only  heat  source  which  influences  the  shape  attained  by  | 

the  sandwich.  j 

• CONCLUSIONS  ? 

1.  Although  plagued  by  reproducibility  problems,  all  possible  combinations  of  \ 

the  four  catalysts  investigated  exhibited  positive  or  possible  synergistic 

effects  on  sample  burn  rate  at  a pressure  of  2000  psia.  j 

2.  The  synergistic  effects  are  weak  and  within  the  range  of  the  reproducibility 
of  the  experiments.  Consequently,  they  are  of  doubtful  practical  utility  with 
the  catalysts  investigated  in  these  experiments. 

3.  The  ferric  oxide-iron  blue  system,  which  was  the  only  system  tested  over  the 
300-2000  psia  pressure  range,  exhibited  weak  posititive  synergistic  effects 
only  at  300  and  2000  psia. 

4.  The  cast  propellant  samples  with  Harshaw  catalyst  Cu  0202  P loaded  in  the 
oxidizer  particles  did  not  burn  as  fast  as  the  binder  loaded  samples  at 
600  and  2000  psia.  However,  difficulties  in  the  experimental  determination 
of  the  actual  catalyst  loading  in  the  oxidizer  particles  render  the  results 
inconclusive, 

5.  Analysis  lias  shown  that  the  faster  are  the  bindor-oxidizer  reactions  the 
more  V-shaped  will  be  a sandwich  when  deflagrating,  reinforcing  prior 
interpretation  of  sandwich  sliapes  and  the  effect  of  catalysis  of  binder- 
oxidizer  reactions  on  surface  shape. 
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APPENDIX  A 

Solid  Propellant  Sandwich 
Deflagration  Analysis  with 
Binder  - Oxidizer  Reactions 


NOMENCLATURE 


constants  in  the  perturbation  solution 

distance  above  solid  surface  at  which  the  mass  fraction  of  H&O^  vanishes 

specific  heat  of  solid  phase 

specific  heat  at  constant  pressure 

deviation  of  c fiom  the  one-dimensional  value 

activation  energy 

dimensionless  temperature,  t/Tq 

deviation  of  g from  the  one- dimensional  value 

specific  enthalpy 

dimensionless  specific  enthalpy  of  formation  h /c  T 

o p o 

specific  reaction  rate  constant 

k pW 

dimensionless  preexponential  factor 
linear  transport  operator 
root  for  exponential  solutions 
mass  flow  rate  per  unit  area 
endothermic  heui,  of  gasification 

heat  of  reaction  in  the  gas  phase  for  the  pure  AP  deflagration 

burn  rate 

temperature 

production  rate,  mass/volune/time 
dimensionless  production  rato 
molecular  weight 
horizontal  coordinate 
vortical  coordinate 
moss  fraction 

mass  fraction  of  IflL  at  solid  surface  for  one  dimensional  case 


Y#  deviation  of  y'  from  one  dimensional  value 

s 

y deviation  of  Y from  Y 

p a + y ' 2 ) 1/2 

a thermal  diffusivity  or  constants  in  perturbation  solution 

0 constants  in  perturbation  solution 

£ dimensionless  activation  energy,  E/RTq 

t\  y - ys 

X thermal  conductivity 

% c X/cX 

P s'  s g 

o density 

t reaction  time 

r 

Subscripts 

cold  solid  values 
reaction  1 or  values  at  T|  = c 
reaction  2 
reaction  3 
reaction  h 
binder 

■» 

HClOj, 

products  of  AP  deflagration 
pi-o ducts  of  reaction  2 
products  of  binder-oxidiser  reaction 
pertaining  to  one  dimensional  AP  deflagration 
gas  phase 

solid  phase  or  at  the  sclid-gas  interface 
in  the  gas  phase  at  the  solid-gas  interface 


0 

1 
2 

3 

k 

B 

P1 

F2 

h 


3 

AP 

g 


Superscripts 

* dimensional  quantity 

' differentiation  with  respect  to  x 

- pertaining  to  one-dimensional  deflagration 
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INTRODUCTION 


In  Ref.  (A-l)  an  analysis  was  conducted  of  the  deflagration  of  a semi- 
infinite  slab  of  ammonium  perchlorate  (AP)  adjacent  to  a semi-infinite  slab 
of  binder.  It  was  assumed  that  the  binder  was  dry  and  passive  so  that  reactions 
between  the  oxidizer  and  binder  were  too  slow  to  be  effective  in  driving  the  de- 
flagration* and  there  were  no  melt  flows  to  complicate  the  binder-oxidizer  inter- 
face region.  Tne  slow  reaction  assumption  required  that  the  AP  self-deflagration 
was  the  only  mechanism  responsible  for  pyrolysing  the  binder.  Accordingly,  the 
analysis  could  only  be  valid  above  the  low  pressure  deflagration  limit  of  the  AP. 

In  the  current  work  the  assumption  of  slow  binder-oxidi2er  reactions  is  relaxed. 
However,  for  analytical  tractability  the  no-melt- flow  assumption  is  retained. 

It  was  found  in  the  previous  treatment  that  regardless  of  binder  type  the 
M‘  would  deflagrate  with  a nearly  horizontal  surface,  for  from  the  binder,  and 
the  surface  would  dip  d‘vr.  sligatly  as  the  binder  was  approached.  There  was  then 
a slope  discontinuity  at  the  binder-oxidizer  interface  and  the  binder,  for  usual 
binder  properties,  would  assume  a nearly  vertical  slops.  By  relaxation  of  the  sere 
binder-oxidizer  reaction  rate  assumption  it  is  desired  t©  find  a modification  to  the 
surface  slspe  that  is  invariably  produced  when  a catalyst  is  added  experimentally 
1 . That  Is,  the  AP  assumes  a stall  acute  angle  with  respect  to  the  nearly 
vertical  binder.  Tne  reactions  between  the  oxidizer  and  binder  drive  the  inter- 
face region  at  a faster  vertical  bum  rate  than  was  possible  when  the  AP  deflag- 
ration above  was  responsible  for  the  overall  sandwich  bum  rate. 

Bgactic-ns  fgnsl  ierr  d 

The  simplest  possible  set  of  gas  phase  reactions  which  it  appears  possible 
to  choose  and  still  show  the  dominant  processes  is  a set  of  reactions  consistent 
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with  Ref.  (A-3)  . This  is 

K 

F1  + F2  ^ P1 
F2 

I'  + B -*  P 

K 


B +Pl^p3 

k4 

P + p -**  p 

2 *1  3 


(1) 


Here  F^  denotes  NH^  and  F0  denotes  HCIO^  presumed  to  be  the  major  reactive 
products  liberated  from  the  AF  surface.  Species  B is  the  binder  pyrolysis  pro- 


duct and  is  some  product  intermediate.  is  the  final  product  of  the  oxidizer 
binder  reaction.  It  is  to  be  noted,  of  course,  that  these  are  a simplistic  set  of 
global  reactions  only  chosen  to  describe  gross  effects. 

In  order  to  apply  fundamental  chemical  kinetics  to  the  above  reactions  it 

wild,  first  be  assumed  that  F^,  F?  and  B have  identical  molecular  weight  W.  Then 

P^  has  molecular  weight  2W,  as  does  P0,  and  has  molecular  weight  3W.  Using 

'A-IA 

fundamental  kinetics  laws  x , the  production  rates  of  the  various  species 
may  be  written 


* 1 <v* 

w = -k^  Yp  ^ k4  YP  XF 

1 fl  f2  * 4 "2  *1 


WP  = kl  YF  YF  ~ k2  YF  YB 
r2  12  l2 


# 

WB  " k2  YFg  YB 


‘ 2 ^3  yb  ypx 


*W>C 


v?l‘  “i  hj  yf2  " k3  h *Pj 


\-K\h  -K  \ \ 


* 


V 3 S YP1  + 2 \ lP2  YF1 


(2) 
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Figure  A-l . Sandwich  Deflagration  Schematic. 


These  expressions  must  be  used  later  in  the  species  conservation  equations. 
Substantial  simplification  will  occur,  however,  if  the  boundary  conditions 
on  the  species  are  first  examined. 


Species  Boundary  Condition 


*)(•  *)(• 

Making  distance  dimensionless  by  0/  /r  and  using  Fick’s  law  of  diffusion, 

s 

as  in  Ref.  (A-l) , the  boundary  condition  at  the  geis-solid  interface  for  any 
species  may  be  written  as 


5^ 

dn  J + 

s 


= -4"  (Y,  - 


.* 

m. 

1 + 
s 


i .* 
s m j 


(3) 


where  d/dn  = (l/z)  ( d/dy)  - (yf/z)  (d/dx) . See  Fig.  (A-l)  for  the  configuration. 

s 

In  order  to  evaluate  the  derivatives  along  the  surface  the  same  approach  as 


used  in  Ref.  (A-l)  will  be  used.  A functional  form  for  Y^  is  chosen  as 


Yi = Yi  + (Yi  ■ Yi  )(y-ys)/c 


(4) 


with  Y.  , y.  , y , and  c as  unknowns  to  be  determined  from  the  conservation 
1s  xl  s 

equations.  All  of  these  unknowns  are  functions  of  x.  Note,  however,  because 

of  the  definition  of  c,  Y„  = 0.  Eq.  (4)  will  allow  evaluation  of  the  dei'ivative 

in  Eq.  (3) . first,  however,  look  at  the  mass  flow  ratio  in  Eq.  (3)  • On  the  binder 

surface  only  binder  can  penetrate  so  the  mass  flow  ratio  is  unity  if  i=B  and  is 

zero  for  all  other  species.  On  the  AP  surface,  consistent  with  Ref.  (A-5) > P^> 

F.  and  F9  are  the  species  for  which  the  surface  is  a source.  Let  G be  the  frac- 

'-ion  of  the  surface  mass  flow  which  i3  P^  and  then  (l-G)/2  will  be  the  mass  flow 

ratio  of  P,  and  F_.  Tills  "dilution  coefficient ",  G,  will  be  assumed  independent 
X 2 

of  x,  in  accord  with  the  independence  on  surface  temperature  shown  in  Ref.  (A-i>) . 
Ml . (3)  cal*:-  for  as  an  example  yields 

’V  £ [V].i[\-^1 

c 2 L S J z U S J 


(5) 


On  the  other  hand,  vixen  evaluating  Eq.  (3)  for  on  the  binder  but  just  at  the 
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interface,  mure  care  must  be  taken.  For  no  mass  sources  or  sinks  the  Y.  and 

h 

Y.  of  Eq.  (4)  must  be  continuous  at  x = 0 and  their  first  derivative  must  be 
s 

continuous.  However  y/  undergoes  a discontinuity  at  x = 0.  therefore,  Eq.  (3) 

s 

just  to  the  right  of  the  interface  is 


1 

“b 


ysB  ,1 


y ' * 


r fi  + ^ Fi  u j.  \ 

L S C S J Z_  2 


(6) 


“B  “ ” B 

where  it  is  understood  that  Eq.  (6)  is  only  good  at  x = 0.  Eq.  (5)  is  a boundary 
condition  that  applies  all  along  - ® £ x ^ 0.  Ihe  boundary  condition  on  B is 


Y Y 
11s 


i r 
z L 


. ys'(\\)  Tm±\ 
C J z 


17) 


but  right  at  x = 0 


Yn  y / 

i - J*  n£  . i (yb  . yb  ) 

Zg  C ZgLSC  1 S 


1 = i-  (yb  -1) 

J Zg  a 

(8) 


Eqs.  (6)  and  (8)  may  now  be  investigated  concerning  orders  of  magnitude. 

A perturbation  solution  will  be  sought  so  that  at  x » 0,  which  is  the  binder- 
oxidizer  interface,  only  a small  perturbation  on  a one-dimensional  AP  deflag- 
ration will  be  taking  place.  That  is,  the  Pg,  P^  and  B mass  fractions  will  be 
small  compared  with  unity  and  the  deviations  of  F^,  Fg  and  P^  mass  fractions 
from  their  one-dimensional  deflagration  values  will  be  small  compared  with  unity. 

Guided  by  Ref.  (A-l)  it  will  be  assumed  that  y*  is  large  compared  with  unity  and 

SB 

y * is  small  but  that  this  product  may  be  of  order  unity.  Derivates  in  x will  bo 
5 

of  order  § compared  with  unity  and  § is  large.  It  is  assumed  that  y*  is  of  the 

SB 

order  of  % or  larger  and  all  perturbation  quantities  are  of  the  order  of  l/§. 

From  Hof.  (A-l) , c is  of  order  l/§.  Using  these  orders  of  magnitude  in  Eq.  (8) , 
tiie  largest  terms , which  are  of  order  §,  are 

(9) 


y'  K (o)  « § 

*B  BS 


Splitting  Yp  and  c into  unperturbed  and  perturbed  parts  and  using  Eq.  (6) , 
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the  boundary  condition  is 


y4  vii0)=  "f^*K7i)-  5yf 


An  equation  similar  to  Eq._(lO)  holds  for  Fg .For  P1  the  equation  is 


y.'  < (o)  = 


SB  P1 


( i + ys'  ys'  ) - § ( 1-2  Yp  ) 

B s 


However,  investigation  of  the  interface  conditions  for  P and  P yield 

* ' (0)  = 0 = Y*  (0) 


Equations  (9)  - (12)  hold  at  x = 0.  Hie  corresponding  boundary  conditions  which 
hold  for  the  unperturbed  and  perturbed  quantities,  along  the  AP  surface,  - ® s x <;  o , 


< 5 *4-)  - S(i.-  s) 

e s 2 


for  the  unperturbed  quantity  and 


C (Y  /c2)  + (5+1/c)^  =0 

s \ 

G (T  /c2)  + (%  + 1/5  ) y_  =0 


(§+i/5)y  . y /s 


(§  + VS)  y„  y /5 


2s  p2. 


C(f  -D/o-^-t 5+1/-c)v  +,  0 (14) 

s s ■L1 

for  the  perturbed  quantities.  In  the  above,  terms  of  order  l/§  compared  to  unity 
liave  boon  dropped. 


Gas  Phase  Solution 


The  conservation  equations  have  been  developed  in  Ref.  (A-l) , so  they  will 
only  be  stated  here.  Using  the  variables  T1  a y - y , x instead  of  y,  x , the 
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energy  equation  is 

^ ( g ) = 

i=d 

i*£+z2  £ -2<  i-  -O-  - si. 

Sx^  3Tf  ST|dx  all  ST) 


(15) 


For  species  i the  conservation  equation  is 


*<*!> 


* * *,.**. 
•i  * 'p  / <r  *g  > 


(16) 


Combination  of  Eqs.  (15)  and  (16)  yields  a simple  homogeneous  equation 

6 

**  ( 8 + r \ Yi  ) = 0 (17) 

i=d  0 


p / 2 

In  Eq.  (15)  the  operator  l contains  z = 1 + y .To  the  accuracy  of  the 

s 

/ 2 

analysis  yg  may  be  neglected  compared  to  unity.  For  the  unperturbed  deflagra- 
tion of  AP,  field  quantities  are  functions  of  11  alone.  The  solution  for  the  AP 
deflagration  has  been  presented  in  Ref.  (A-l) , except  for  the  modification  due  to 
the  dilition  coefficient  which  is  included  in  this  work.  The  solution  is 


yf  < = + 5 ) = ? ( 1-G  ) * Q 

1 0 IT 

s 1 


Q s \ 4 S I )2  e^P  { - ei  /Cia  + (i,.  - Sg  ) ? ] }d(5-) 

1_  0 c 1 c c 


h ' ®s  ~ 2 % YI\ 


% 6 *F. - 


1.  h 


(18) 


In  this  work  the  thermodynamic  properties  of  F^  and  Fg  are  assumed  identical  to 
each  other,  for  simplicity.  All  numerical  computations  have  been  performed  for 
?=  9*0  , =4.022,  gs  = 2.93,  G ~ .75>  c = .119  and  Yj.(,  = .0634.  This  corresponds 

to  burning  at  800psia,  but  the  dimensionless  results  are  squite  insensitive  to 
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pressure,  as  explained  in  Ref.  (A-l) . 

The  solution  of  Eqs.  (14)  and  (15)  is  by  an  integral  technique.  Using 
Eq.  (4)  as  a guessed  form  for  the  mass  fractions  and  g - g a(g  -g  ) T/c  for 
the  temperature  profile,  Eqs.  (16)  and  (17)  are  integrated  in  T|  from  0 to  c. 
This  yields  ordinary  differential  equations  in  x.  However,  appearing  in  these 
equations  from  integration  of  the  second  term  in  l from  Eq.  (15)  is  the  T1  - 
derivative  of  the  appropriate  quantity  evaluated  at  T]  » c.  This  is  an  unknown. 
To  eliminate  this  unknown,  as  discussed  in  Ref.  (A-l) , the  parabolic  form  of 
Eqs.  (14)  and  (15)  are  used  at  T]  = c.  That  is,  in  Eqs.  (l4) . 


fh  \ - 5 

^2  ) T)=c  w ST\=c 


is  used  to  eliminate  dY^/dl]  ) ^__c  in  favor  of  the  x - derivative  and  reaction 
rate.  A similar  operation  holds  for  Eq.  (17) . 


The  next  step  in  the  solution  is  to  let  each  unknown  be  represented  by  its 
unperturbed,  one-dimensional  value  plus  a perturbation  component.  The  equations 
are  then  split  into  unperturbed  equations  and  perturbation  equations  and  the 
perturbation  equations  are  then  linearized.  The  unperturbed  equations  have  Eqs. 
(18)  as  the  solution,  after  using  Eq.  (13) . 


It  is  first  instructive  to  view  the  solution  for  the  binder  mass  fraction. 
The  differential  equation  is 

yBs  (c/2)  + yBx  (c/2  + 1/s) 

+ 'VS  )<  5 + vs ) - £2  { ''Bs  h **\  \ } 

s 1 

+ 2 *3 1 { yBs  \ + ^ \ } + i k3  v\  ^3/sl  (19) 


where  the  reaction  rate  integrals  are  evaluated  from  the  unperturbed  solution 
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2S  - J YFs  (1-H/5)(1-Vc)  e'Vg  . 


dn 


:2X  = J c yf  ( i - yi ) e_Vg  dn 

o 

3S=]  \0l)  ( 1 - Vc  ) e"e3/g  dTl 


I31  = j ° \ e“£3/g  <B\ 

o 


Notice  that  all  I’s  are  positive.  Now  viewing  the  hinder  boundary  condition  in 
Eqs.  (l4)  j yfi  06  Consequently , Eq.  (19)  is  a homogeneous,  linear  second 


order  differential  equation  with  the  solution 


\ - \ • v 


Making  the  required  substitution  bl,  , . . 

a satisfies 


4 ■ip^jjV  v).*e/S 

(f  + 1/S  ) (c/2  + l/§  ) + 1/2 


(20) 


where 


t'hh  * 


H 


v-*8  Vs*3 


2 

By  inspection, therefore,  is  positive  and  ra^  has  two  real,  equal  roots  of 
opposite  signs,  the  positive  sign  is  accepted  so  thatyB  will  decay  to  zero 


at  x - -«•.  The  constant  of  integration,  Ag  is  determined  by  application 
of  Eq.  (9)  S 


\ =t/(v;B> 


The  binder  solution  is  therefore  completed  without  knowledge  of  the  behavior 
of  the  other  mass  fractions.  Either  reaction  2 or  reaction  3 will  increase 
[see  Eq.  (20) ] as  the  reaction  rates  increase.  This  has  the  effect  of  lowering 


'3  n 


1 1 


1.-  hVaer-V.f  Ti  V^gl 
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the  amount  of  binder  present  in  the  gas  phase  [Eq.  (21) ] and  of  increasing  the 
decay  rate  of  binder  as  x proceeds  away  from  zero.  These  are  obvious  physical 
facts  since  reactions  2 and  3 are  disappearance  reactions  for  the  binder  gases. 
Also  note  from  Eq.  (20)  that  ra^  is  indeed  order  §,  if  £ is  less  than  or  equal 
to  order  unity.  As  a consequence  Yl  is  of  the  order  of  §/_,  as  asserted  above. 

D 13 

For  £ = y=  0,  ra  = 6.73*  Note  furthermore  that  m^  is  independent  of  y^,  under 

the  restriction  already  imposed  that  "y / is  of  order  l/§. 

s 

All  other  mass  fractions  do  not  behave  in  an  autonomous  manner,  however. 

The  differential  equation  derived  from  continuity  of  the  Pg  species,  for  example, 
yields 


\ (5/2)  + \ (5/2  +1?) 

s 1 


+ ( 


y )(  5+1/e) 

P2, 


1 

/ 


+ k, 


{ 


\ + V \ * 

s 1 


(22) 


where,  as  before, 

;>C 


I.  = ! 

l » 

s o 


X 00(1-1/5)  e’V* 


dH 


\ = J Yf  (D  e"V*  d'Tl 

1 o 

Considering  that  V * V from  Eqs.  (Ik) , Eq.  (22)  is  a linear,  second  order 

p2  p2. 

s 1 

equation  in,  say,  Pg  , but  it  is  inhomogeneous  because  of  the  appearance  of  the 

s 

binder  mass  fractions.  It  follows,  therefore  that  Y will  have  a solution  of 


the  form 


- ^ 


2s  e 


P2X  + raBx 

Pg  e 
s 


It  may  be  verified  by  substitution  that  the  root  to  the  homogeneous  equation, 
mp  , increases  with  an  increase  in  rate  of  the  destruction  reaction,  reaction  4. 


Application  of  Eq.  (12)  yields 


= - a 


m 


Continuing  for  F^,  the  solution  would  look  like 


i/  . m 

%.  - V ep3  +a 


, -V*»p  .V 

*3  *3 

Js  ^s 


since  the  production  of  depends  upon  both  B and  Pg.  However,  since  P^  is 

never  destroyed  m is  determined  hy  Eq.  (20) , setting  £ = y =0. 

P3 

Here  m = m = 6.73  an£  all  decay  scales  other  than  a fixed  m are  set  by  m 
p3  p2 

and  nig. 

The  situation  with  F^Fg  and  P^  are  more  complex,  however,  since  they 
contain  unperturbed  as  well  as  perturbed  components.  Only  two  of  them  need 
be  considered,  however,  since 

y +y+y  + V + V = o 

F1  F2  h *2  P3 

Here  will  be  eliminated  and  only  F^  and  Fg  are  considered.  Furthermore, 
consider  now  the  case  where  kg  and  k^  are  zero  (only  a diffusion  flame  be- 
tween the  binder  and  AP  deflagration  products  exists) . In  this  cave  the  eq- 
uations for  Fj^  and  Fg  are  the  same , as  are  the  boundary  conditions . 

Consequently  ^ , as  is  obvious  since  there  are  no  chemical  reactions 

1 £ 2 

that  prefer  one  over  the  other. 


The  differential  equation  for  conservation  of  F^  (or  Fg)  is 

C#/(Yf/2)  + VT"(l/2)  + V?  /l  - CYF  /c2 
s Xs  1s  S 

+ y"  Yp  + § y,  = (bqjb yf  ) ^ + (dQ/dc)  c 

S 1S  1s  1s 

+ (3Q/5y1)  G1  + (SQ/9gs)  Gg  (23) 

where  the  derivatives  of  Q may  be  evaluated  from  Eq.  (18) . The  unknowns  are  now 

C,  y F , y|  and  Gg.  The  boundary  condition  of  Eqs.  (l4)  provides  an  algebraic 
Xs 
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relation  between  C and  Yp,  > A second  algebraic  condition  comes  from  the  assump- 
tion of  an  equilibrium  gas-solid  interface  , whereby 


F1  = ~Fe  "eB'/Ss  , \ _ °,  eB 


G e. 


(24) 


Ss2 


A fourth  equation  in  the  five  unknowns  is  obtained  by  integrating  Eq.  (15)  and 
applying  the  perturbation  procedure 


(0i"  + Gs"  )(S/2)  + (os  - Gt)(  5 +1/5 ) 

+Gi"/5+y~  ^ - a ){ (v,"  +n"Hs/2) 

Z o P1  L X1  XS 

f1,f?jp3,b  0 

+ * ( *i  - ^ X E + 1/5  ) } = 0 (25) 

1 S 1 

Equation  (25)  is  nut  a homogeneous  equation  in  the  five  unknowns  because  under- 
neath the  summation  sign  sit  ^ and  y^,  whose  solutions  are  now  known. 

In  the  current  case  the  operation  on  F in  Eq.  (25)  is  the  homogeneous  operator 

mPox  nix  ni.x 

which  wipes  out  the  e = e solutions , One  is  therefore  left  with  only  e a 

type  of  behavior.  The  only  way  for  the  four  equations  in  five  unknowns  to  have 
a solution  and  still  satisfy  the  boundary  condition  of  Eq.  (10)  is  to  let 


G-,  = A e 

1 Gi 


mux 


> C = Afi  e 


infix 


msx 

ffiQX 


•°s  = V 


Y Y 

\ = % = 
s 5 


max 


and  use  the  boundary  condition  as  an  additional  equation.  It  is  to  be  noted  that 

this  procedure  also  holds  if  ••=  0,  the  case  of  a purely  inert  binder.  This  latter 

case  is  the  one  treated  by  Ref.  (A-l) , but  the  solution  there  did  not  involve  Eq. 

(9) . It  is  now  believed  that  Ref.  (A-l)  is  in  error  since  Eq.  (9)  must  be  satisfied 

for  consistency.  In  any  event,  following  the  solution  procedure,  eliminating  C and 

G ani  recalling  that  in  Eq.  (10)  y*  (0)  = y;  + A(  , there  result  three  equations 
* s s y 
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in  three  unknowns  of  the  form 


*12  *y 


F 


*21  \ + *22  *Yf 


+ A 


23  *Y  ‘ 


0 


*32  *Yp  + *33  V 


(26) 


It  may  be  verified  in  Eq.  (26)  that  A^<0  and  for  exothermic  reaction  between  B 
and  P^j  the  heat  of  formation  difference  in  the  first  of  Eq.  (26)  is  negative. 
Since  £ and  y are  positive,  the  first  of  Eq.  (26)  clearly  shows  that  the  re- 
action tends  to  elevate  g^. 


RESULTS 


Consider  first  the  solution  for  no  reactions,  £ = T = 0.  For  various  y / 

/ Sb 

and  y*  values  the  solutions  to  Eq.  (26)  are  shown  on  Fig.  A-2  as  A vs  y*  (o)= 
-/  / , 61  S 
y + Av- . The  values  for  y*  increase  as  one  moves  to  the  right  along  any  curve 
s 1 w f SB 

of  fixed  y . It  is  seen  first  of  all  that  there  are  no  solutions  for  which 
s 

yl  (0)  >0.  Furthermore , in  regimes  where  - g^  are  leas  than  zero,  which  is 
the  only  physically  reasonable  regime,  since  the  binder  is  a heat  sink,  y*  <y#. 
This  means  tiiat  the  AP  surface  turns  down  to  attempt  to  intersect  the  binder  at 
an  acute  angle,  rather  than  an  obtuse  angle.  This  is  demanded  by  the  gas  phase 
solution  alone  and  does  not  depend  upon  binder  or  AP  properties  in  the  solid 
phase.  This  character  of  the  solution  is  somewhat  counter  to  the  author's  in- 
tuition. It  was  expected  that  the  AP  would  curve  upward  as  the  binder  is  ap- 
proached . 

For  the  cose  with  reaction  consider  a high  activation  energy  case  so  that 
Y is  dominated  by  the  last  terra  and  £ < < y . As  an  example  takey  » 0.15  and 
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figure  A-2.  BsweSwe  at  the  Bags  of  ti:c  Ar 

" »_f* f ' “2  "f  “» I'-Tace  ^^t 

tiie  ainier-Oxidisei'  Interface. 


3^ 


S.  - fi  = -10.  On  Pig.  A- 2 the  dashed  line  shows  the  solution  for  y„  =0.1. 
P3  PX  B 3 

The0temperature  is,  of  course,  elevated  as  compared  with  the  no  reaction  case. 


Here  there  is  no  good  physical  argument  to  reject  solutions  for  A >0.  However, 

gl 

only  one  solution  curve  is  shown  and  for  yg  =»  0 .1 , A^  < 0 . 


In  order  to  show  which  binder  properties  belong  to  each  point  on  the  curves 
of  Fig.  A-2  one  must  consider  first  the  solution  for  the  solid  phase  heat  transfer. 
This  solution  was  obtained  in  Ref.  (A— 1) . By  inspection  of  that  solution,  invoking 
the  continuity  of  the  heat  transfer  vectors  in  the  solid  and  gas  phases  at  the  in- 
tersection of  the  binder  and  oxidizer,  and  applying  the  order  of  magnitude  argu- 
ments, the  solution  for  the  binder  heat  of  gasification  is 


% = “a  As  ys  * [ 1 * <.  ys  <0)  +1-«J 

S Zj 


3 


(27) 


§ c 


It  turns  out,  therefore,  tiiat  the  result  is  independent  of  the  solid  phase  solution, 
mis  is  a remarkable  property  of  the  problem.  The  final  link  to  the  binder  is  the 
pyrolysis  law.  As  stated  in  Ref.  (A-l)  this  is 


1 


(23) 


Given  the  activation  energy,  b„  is  directly  determined  by  y*  , since  r <=  r, 

8 

consistent  with  the  orders  of  magnitudes  considered  here.  Then  • (2?)  deterstinea 
tiie  allowable  for  any  gives  y*,  which  is  the  eigenvalue  of  the  problem. 

the  results  of  this  calculation  are  plotted  for  the  no  reaction  case  in 
rig.  A-3.  Given  a set  of  binder  properties  ieterainea  the  eigenvalue  y*  on  this 

9 

plot.  She  disappointing  result  is  that  within  the  solution  constraint  that  y*  is 
araa.il  there  are  no  solutions  for  binder  properties  of  interest.  Shown,  fsr  example, 
are  the  CUB  and  HtSfi  points  for  properties  taken  froa  Kef.  \A-6) . There  siuply 
are  no  solutions  to  the  problem  for  which  < 0 at  large  positive  q^.  Vhat 
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® HTPB  at  800  psia 
&CTPB  at  800  psia 


tion  C i 


must  occur  for  binders  of  interest  is  one  must  consider  the  fully  nonlinear  set 
of  equations.  Nevertheless,  this  solution  may  be  used  as  a guide  to  investigate 
the  qualitative  effect  of  parameter  changes. 

Within  the  regime  of  allowable  solutions,  which  basically  limits  the  solution 

to  small,  positive  or  negative  q^,  and  y / >0,  in  increase  in  pyrolysis  rate  at  fix- 

ed  q_  increases  y/.  An  increase  in  heat  of  gasification  at  fixed  pyrolysis  rate 
ij  s 

(y ' fixed)  decreases  the  slope  of  the  AP  far  from  the  interface  (y7  decreases) . 

SB  S 


For  the  case  with  reaction,  considered  above,  the  solution  for  y7  = 0.1  is 

s 

shown  on  Fig.  A-3.  Exothermic  reaction  shifts  all  of  the  curves  to  lower  qn  for 

fixed  y 7 . Therefore,  for  a fixed  set  of  binder  properties  the  exothermic  reactions 
®B 

decrease  yg.  That  is,  the  solution  says  that  reactions  tend  to  promote  an  acute 
angle  between  the  oxidizer  and  binder,  as  is  experimentally  observed. 

Tracing  through  the  nondimens ionali zation  procedure  and  for  the  numbers  used 

here  a characteristic  reaction  time  is  given  1/r  = e J^p/w  , 

^ 3 

and  t alO  sec  at  800  psia  for  y„  55  0.1  (which  depends  on y.  through  Eq.  (9) . 

r SB 

From  Ref . (A-7)  a typical  reaction  time  under  these  conditions  may  be  computed  to 

be  10  ^sec.  Therefore,  typical  reaction  times  are  well  within  the  range  tnat  a 

significant  impact  would  be  made  on  the  sandwich  shape  by  the  binder-oxidisor 

reactions.  Furthermore,  for  completeness  reactions  2 and  4 sliould  be  considered. 

However,  because  of  the  limited  applicability  of  this  solution  together  with  the 

qualitative  information  already  gleaned  from  the  solution,  no  further  computations 

will  be  made. 

CONCLUSIONS 

1.  A linearised  solution  using  an  integral  technique  has  been  found  to  the  problem 
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of  a deflagrating  semi-infinite  slab  of  AP  adjacent  to  a semi-infinite  slab 
of  binder  which  is  capable  of  reaction  with  the  AP  in  the  gas  phase.  Hie 
solution  was  found  to  only  be  valid  for  binders  which  are  difficult  to  gas- 
ify (low  pyrolysis  rates)  and  with  low  or  negative  (exothermic)  heats  of 
gasification.  A further  restriction  to  the  solution  is  that  there  are  no 
melt  flows. 

2.  For  typical  reaction  times  the  binder-oxidizer  reactions  are  fast  enough  to 
substantially  affect  the  deflagration  rate.  The  effect  is  to  make  the  slab  of 
AP  intersect  the  binder  at  an  acute  angle,  giving  the  appearance  that  the  in- 
terface region  is  driving  the  deflagration  rate,  which  it  is. 

3.  The  solution  predicts  that  as  one  approaches  the  binder  along  the  AP  surface 
the  AP  surface  becomes  convex  from  above,  which  is  somewhat  counter  to  physi- 
cal intuition. 
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